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Abstract 

We propose a new method for the design of adaptation algorithms that guarantees a certain prescribed 
level of performance and applicable to systems with nonconvex parameterization. The main idea behind 
the method is two-fold. First, we augment the tuning error function and design the adaptation scheme 
in the form of ordinary differential equations. The resulting augmentation is allowed to depend on state 
derivatives. Second, we find a suitable realization of the designed adaptation scheme in an algebraic-integral 
form. Due to their explicit dependence on the state of the original system, such adaptation schemes are 
referred to as adaptive algorithms in finite form, in contrast to (conventional) algorithms in differential form. 
Sufficient conditions for the existence of finite form realizations are proposed. It is shown that our method 
to design algorithms in finite form is applicable to a broad class of nonlinear systems including systems with 
nonconvex parameterization and low-triangular systems. 
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1 Introduction 

Significant progress in adaptive control theory has been made in the areas of linear and nonlinear systems p^ ll.'il 
1201 EHl I plants with relative degree greater than one |1(J[I19|[T^ . and systems with nonconvex parameterization 
|17[ I15|. However, there is still room for further developments, as there are important unresolved problems 
regarding the issue of performance, especially in the presence of nonconvex parameterization. 

As expressed by asymptotic stability of adaptive systems (201 j robustness, and good transient behavior |14j. 
suitable performance can be proven under the requirement of persistent excitation. As it is generally observed 
in practice, insufficient excitation results in absence of asymptotic stability and, as its consequence, in poor pa- 
rameter convergence, sensitivity to small disturbances and poor transient performance. A suitable performance 
criterion is needed to assure the efficiency and quality of the system. As substitutes for performance criteria, 
most of the available results in direct adaptive control without restrictive persistent excitation requirements 
limit themselves to L2 and Lqo^ norm bounds of the tracking errors. For more sophisticated performance mea- 
sures like the LQ criterion, some results are available [71|B]. These results, however, deal either with too narrow 
a class of uncertain systems or present only a comparison between adaptive and robust backstepping [S] 
without suggesting new adaptation schemes. On the other hand, when improvement heuristics are suggested 
like in |21| , no exact performance criterion is provided that can explicitly be computed a-priori, except probably 
the bounds on L2 and Loo norms for the tracking errors. 
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Another unresolved issue in conventional adaptive control theory is nonconvex parameterization of the plant 
model. The available approaches encourage to compensate for the nonlinearity by using an additional damping 
term These techniques guarantee existence of a solution to the control problem, but are limited in 

terms of their practical value because of their high-gain nature. Moreover, none of these approaches provides a 
performance measure beyond Lac or L2 norm bounds with respect to the state or tracking error. 

An important impediment to solving these two major problems at once, we believe, is the lack of sufficient 
information in the conventional adaptive schemes. One way to provide the algorithms with extra information is 
to augment the tuning errors. This idea is inherent to both Morse's adaptive controllers ^Hl and those based on 
Kreisselmeier's observers jl(J) when dealing with plants with relative degree greater than one. These augmented 
errors then are used in conventional gradient schemes 

6» = -rV'(x,t)^(x,0,i), (1) 

where x G i?" is a state (or output) vector, e i?'' is a vector of controller parameters, -4(x, 6) is an operator 
that depends on the particular problem at hand, V'lx, t) is the error function and gain F > 0. Improved 
performance of these augmented controllers for Morse's high-order tuners is reported in In these papers 

the plants are assumed to be linear, but no additional criteria are provided except for L2 and Loo norm bounds 
with respect to the tracking errors. 

The above-mentioned limitations of performance and applicability motivate us to search for a new augmen- 
tation that uses additional information about the system dynamics, such as, e.g., state derivatives. As a result, 
new properties in the system can be created. On the other hand, we wish to find physically realizable algorithms 
that do not require measurements of any unknown signals, derivatives, or parameters. In order to meet these 
seemingly contradictory requirements we propose to extend conventional algorithms as follows: 

0(x(t),t) = 0p(x,i) -K07(t); 6»/ =^2(x,0,t), 0p(x,i) = A(x,0- (2) 

Notice that if functions 6{x.,t) are written in differential form (Q, they may depend on unknown parameters 
and unmeasured signals, e.g, state derivatives. Thus, the equivalent description of adaptive algorithm |2Jl in 
differential form may produce an augmentation that is in fact derivative-dependent, thereby providing the 
algorithm with more information about the plant uncertainties. These observations lead to quite unexpected 
consequences. Instead of restricting the design procedure to the algorithms given in differential form it 
becomes possible to design adaptation algorithms in two steps. First, search for the desired augmentation to 
obtain the required adaptive control properties. At this stage it does not matter whether the augmentation is 
uncertainty-dependent or not. Once a suitable tuning error is chosen, second: find a realization of the algorithm 
in the form of integral-algebraic equation ^ . Such a realization will be termed algorithm in finite form. 

According to our knowledge, algorithms ^ have been introduced for direct adaptive control of nonlinear 
systems in 1986 in ^ and then were reintroduced later in |24l IT] . Their distinctive performance properties and 
extended applicability, however, were not appreciated at the time. The efficiency of the proposed algorithms in 
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was limited by restrictive pseudo-gradient assumptions on 0p(x, t) (for the details see Nevertheless, 
it has been reported recently that algorithms (0) may be able to deal with nonconvex parameterization (see 
for example [301, Lemma 1, p. 558; I^EHI) and guarantee improved transient performance |27||2H|- Some 
preliminary analysis of the distinctive properties of algorithms ||2J) is available, for example, in |2Hll2ni- Ad- 
ditional support and motivation for algorithms (0) can be found in 1 , where the authors introduced their 
adaptation schemes from immersion and invariance principles. However, the main problem, with the current 
study of algorithms l(2Jl is that there is no systematic method that allows us to design these algorithms with 
guaranteed improvements in performance and, at the same time, achieve applicability to systems with nonlinear 
parameterization for a sufhciently broad class of nonlinear dynamical systems. 

In our present work we suggest a new method to design adaptive algorithms in finite form (j^J that guarantee 
improved performance and in addition are applicable to a class of nonlinear ly parameterized plants. The method 
is systematic and is based on two fundamental ideas in adaptive control theory: augmentation of the error and 
embedding the original system dynamics into one of a higher order. These ideas are embodied in two independent 
stages of the design. The first stage is augmentation of the tuning error for algorithms in the conventional 
differential form Q in order to ensure improved performance and extended applicability of these algorithms. 
The resulting augmentation may not necessarily be independent on the uncertainties or time-derivatives of the 
state vector. The resulting augmentation, however, should guarantee certain desired properties of the system. 

Based on the augmentation (possibly, derivative-dependent) obtained in the first stage, the second stage 
of the design method should be to find functions ^i(x, <), A2{x,d,t) which guarantee that algorithms (|2J) 
realize the desired adaptation scheme. We show that this problem may require finding a solution of a system 
of partial differential equations. It is well-known that such a solution may not exist in general. To avoid this 
problem, we consider several special cases of plant models, with their structures satisfying sufficient conditions 
for the existence of a solution. As soon as these basic structures are found, we embed the original system 
into a system of higher order for which the solution is known to exist. The embedding is to be made in such 
a way that the extended system belongs to one of already established basic classes that guarantee existence 
of the solution to the realization problem. With embedding we shall be able to obtain adaptation schemes 
that guarantee not only square integrability of the error but also integrability of its first derivatives as well as 
square integrability of control efforts injected into the system due to the parametric uncertainties. In addition, 
we provide the conditions for which the decrease of the parametric uncertainties and exponential convergence 
into a neighborhood of the target manifold are guaranteed without the restrictive assumption of persistent 
excitation. Last but not the least, our new adaptive schemes can also be applied to systems with nonlinear 
parameterization. 

The layout of the paper is as follows. In Section 2 we specify the class of nonlinear dynamical systems under 
consideration and select the desired augmentation. Section 3 contains the main results of the paper. In Section 
4 we present an example of the design and results of computer simulations. Section 5 concludes the paper. 
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2 Problem Formulation and Preliminary Results 

Let the plant mathematical model be given as follows: 

Xi = /i(x) +5i(x)M, i = l,...,m 

Xj = /j(x) + j/j_,„(x, 0) +.gj(x)u, j = TO+l,...,n, (3) 

where x £ i?" is a state vector, fi, gi : i?" i?, G C^, G ile C i?'* is a vector of miknown parameters, 

Vi : i?" X i?"* — > i?, j/j e C^, u is a control input. Let us define functions f(-), g(-), v{-, •): 

f (x) = (/i(x), . . . , /„(x))^, g(x) = (gi(x), . . . ,5„(x))^, z.(x, 6) = (i.i(x, 0), . . . , J/„_„(x, 6)^ 

It will be useful sometimes to think of state vector x G £ C i?" as x = Xi ©X2, Xi = (xi, a;2, . . . , Xm)"^ , ^2 = 
{xm+i,Xjn+2; ■ ■ ■ , Xn)^ , where symbol © denotes concatenation of two vectors xi e £1 C i?™, X2 G £2 ^ i?"^™ 
and C, Ci, C2 are linear spaces. The time-derivative of xi is independent on 0, whereas the time-derivative 
of vector X2 depends on unknown parameters explicitly. Therefore we refer to the spaces Ci and C2 as 
uncertainty-independent and uncertainty- dependent partitions of system (PJ, respectively. To denote the right- 
hand sides of the partitioned system, we use the following notations: fi ~ (/i, . . . , fm)^ , ^2 = {fm+i, ■ ■ ■ , fn)^ , 
gi — {gi, . . . ,gm)^ 1 g2 — (Sm+ii ■ ■ ■ T9nf^ ■ Hcncc, the partitioned system can be written as follows: 

ii = fi(x) +gi(x)u, ^2 = f2(x) + z/(x,6») +g2(x)u. (4) 

In analogy with the definition of independence of a function with respect to the components Xi of its argument 
x, we would like to define a notion of independence of the function with respect to the partition. Let x G 
C — Ci ® £2, and let fmiction a;(x) : i?" i?" be differentiable for any x G i?". Function a;(x) is said to be 
independent on partition C2 iff 9a; (xi ©X2)/9x2 = 0. We would also like to extend the standard definition of 
the Lie derivatives to the partitioned system. Given the following partition x = xi © X2 ® ■ ■ • © x,. , we denote by 
symbol Lf.ilj{x.) the following derivatives {dip{x)/dxi) fi(x) = if.?/;(x), where f,;(x) stands for the corresponding 
vector- function in f(x) = fi(x) ©f2(x) © ■ • • ©fr(x). If function i/'(x) = (■!/;i(x), . . . , ^„(x))'^ is a vector- function 
then symbol Lf.ip{x) denotes the following vector: Lf.ip{x.) — (Lf.?/'i(x), . . . , Lf;-0„(x))"^. 

As in |3U[ I25| . we define the control goal as reaching asymptotically a target manifold. We assume that 
the target manifold can be given by the following equality ip{x,t) — 0, where ip : i?" x R ^ R, ip{x,t) ^ C-^ . 
Additional restrictions on the function ip{'x.,t) are formulated in Assumptions 

Assumption 1 (Boundedness of the Solutions) Function ip{x,t) is such that for any 6 > there exists a 
function e : R+ —^ R+ such that \ijj{x,t)\ < (5 ||x|j < e{d) along system ^ solutions. 

Assumption^simply states that any trajectory of system ^ belonging to a neighborhood of the target manifold 
'0(x, t) is bounded. Clearly, most of the common goal criteria used in adaptive control satisfy this property, for 
example, positive-definite functions V'(x) for nonlinear systems and quadratic forms for linear ones. In general. 
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however, in order to show the boundedness of x it is not necessary for the function 'ip{'x.) to be positive definite. 
As an illustration, consider the following system: 

ii = fiixi, . . . ,Xi) + Xi+i, i ^ {1, . . . ,n - 1} 

in = /„(a;i, . . . ,a;„) + z^(x, 0) + w. (5) 
Let V'lx) = Xn— p{xi, . . . , Xn-i) + /n-i(a;i, ■ ■ ■ , Xn-i), fi{-),p{-) E C° and furthermore, let the system 

Xi = ft{xi, . . . ,Xi) + Xt+i, i ^ {1, . . . ,n -2} 

Xn-l = p{xi, . . . ,Xn^l) + V (6) 

state be bounded for any v G L^c (i-C- system has the bounded input - bounded state property). Then for 
the system of equations ((S)), it is sufficient that system ^ is input-to-state stable with respect to input v to 
satisfy AssumptionQf with V'(x) — Xn — p{xi, . . . ^Xn-i) + /„_i(a;i, . . . ,a;„_i). 

Assumption 2 (Regularity) For any x € i?" and t > functions -(/'(x, t) and g(x) satisfy the following 
inequality |Lg'0(x, t))| > > 0. 

Assumption El ensures the existence of feedback that transforms the original system into that of the error 
model with respect to the variable ip{x,t). Let Assumption |5] hold; consider 

^ = if^(x, t) + i,(,,e)V'(x, t) + (LgV(x, t)) u + 5^(x, t)ldt. (7) 

Because of Assumption |21 there exists the control input 

«(x, e, t) = (ig(x)^(x, t))-^ (-(p(V') - LfV'(x, t) - i,(,,e)^(x, t) - 9V(x, t)/dt) , (8) 

where G fig C i?'^ - is a vector of controller parameters that transforms {Tj) into 

= -(^(tA) +2(x,6>,i) - z(x,0,t) (9) 

where z(x,9,t) = Li^(x,e)V'(x, i). Let the closed loop system satisfy some additional requirements: 

Assumption 3 (Certainty Equivalence) For any E i}g there exists 6 £ i7g C R , such that for all 
X E R" , t E i?+ the following equivalence holds 

+ i.(x, e) + g(x)«(x, e*,t))] + + ^(^) ^ ^(x, 0, e*,t) + ^(v) = o. (lo) 

ox ot 
It is clear that if Assumption [3 holds then Assumption |21 is automatically satisfied. According to Assumptions 
|21and|3| it follows that z(x, 6 ,t) ^ z(x, 6, t) for any x E i?" and time t > 0. 

Assumption 4 (Stability of the Target Dynamics) Function ip{tp) in flU^) satisfies 

(fiiij) E C", > VV- 0, lim / ip{^)d<, = oo. (11) 



^For plants with uncertainties in functions fi{-), the time-dependent target manifold il'(x,t) = may be required for control 
design. This will be clarified later in the proof of Theorem |1] and in the examples section. 



Assumption 5 (Monotonicity and Linear Growth Rate in Parameters) There exists function a{x,t) : 
R"- X such that (z(x,0,i) - z(x, 0*, t))(a(x, t)^(0 - 6*)) > V z{x,e*,t) ^ z{x,e,t). Furthermore, 

the following holds: \z{x, e,t) - z{x,e* ,t)\ < D\ci{yi,t)'^ {9 - e*)\, D > 0. 

Assumption |21(certainty equivalence or matching condition) simply states that for every unknown 6* G fig 
there exists a vector of controller parameters 6 (9*) G fig such that the system dynamics with this control 
function satisfies the following equation = —ip^ip). Assumption 0] specifies the properties of function (p{ip), 
thus stipulating asymptotic stability of manifold V'lx, t) — for 9^9 and ensuring unbounded growth of 
integral ip{<;)d(^ as —^ oo. Assumption |5l is given to specify the admissible nonlinear parameterization of 
the controller. For linearly parameterized plants this assumption is automatically satisfied. Sometimes we will 
further restrict the nonlinear in parameter functions by: 

Assumption 6 There exists a positive constant Di > such that for anyx., 9,9 ,t>0 the following inequality 
holds: \z{x,9,t) - z{x,9*,t)\ > Di\a{y.,tf {9 - 9*)\. 

Throughout the paper we will assume that functions q:(x, t) and u(x, 9, t) are both bounded in t. For the sake 
of convenience and if not stated overwise we will also assume that functions q:(x, t), '(/;(x, i) are differentiable 
as many times as necessary if differentiation is required to design the algorithm. In addition, we will use 
the term "smooth functions" to denote those functions that belong to C°° . Though it is not necessary at all 
for us to require existence of infinitely many derivatives of the functions that we refer to as smooth in the 
paper, this notational agreement will free the presentation of numerous insignificant details in the formulations. 
Furthermore, along with already defined L2 and L^o norms we will use the following notation: \\9 — 9{t)\\^_i — 

[9 -9{t)rT-\9 -9{t))- 

As mentioned in the introduction, we propose to design the adaptive algorithms in two steps: 1) search for 
the suitable augmentation ensuring the desired properties of the control system, and 2) find the appropriate 
realization of this algorithm in finite form. Therefore we start with the choice of tuning errors ip{'x., t) and 
operators ^(x, t) for the class of algorithms given by formula 9 — r-i/;(x, f)^(x, 9). As a candidate for the 
augmented error ip{K,t) we select the following i) = "0 + tp{x,t). It has been shown in |25[ I3U| that the 
algorithm 

9 = r{ip{iP) + V')a(x, i), r > (12) 
with control ^ guarantee that ?/'(x(t), i) — > as i — + 00 for the closed loop system 

Xl = fi(x) + gi(x)u, X2 = f2(x) +I/(x, 0) +g2(x)'U 

9 = r{<f{i;)+i;)a{x,t), r >0 (13) 

under Assumptions ^ El ~ El In addition, it is possible to show that system ^I'A^ has better performance than 
that of the known schemes. This follows from the next theorem (see also Proposition 1 below): 

Theorem 1 Let system be given and Assumptions\^\^hold. Then for system the following hold: 
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P2) \\6 — ^(i)|jp_i is non-increasing; 
P3) z((x, 6», t) - z(x, e{t), t)) e L2. 



Furthermore, 



Mml < 2Q{i^) + \\eio)~e*\\f,^^y,, \ml<2Qii^) + \\m'0*\\f2Dr)-^ 

llV'lloo < A(Q(V') + ||0(O)-r|| 



(14) 



where Q{')jj) = /(- 



■V(x(0),0) 



ip{(,)d<; and A{d) = max|^|{|?/'| | (y9(<;)(i<; = d}. 



// ^ssMmption is satisfied and function z(x, 0,<) is locally bounded with respect to x, and uniformly 
bounded with respect to t, then 

P4) trajectories of the system are bounded and ?/'(x(t)) ^0 as t ^ 00; 

If in addition functions ip, z(x, 9, t) G ; derivative 9z(x, 0, t)/dt is uniformly bounded in t; function a(x, t) 
is locally bounded with respect to x and uniformly bounded with respect to t, then 
P5) ijj ^ as t ^ 00; z((x,0*,t) - z(x,0(i),<)) -^0 as t^ 00. 

The formal proof of the theorem is given in Appendix 2"^. 

Notice that the function ip{')p) is nonhnear and that its shape influences the L2 and Lao norm bounds for 
^p and tp. Because of this, according to p4|l it is possible to improve the performance of the system with 
respect to L2 and ioo bounds by varying the function ^p{'>jj). The bounds obtained for the Loo norms may be 
improved further for the case when the fmiction ip{ijj) is linear in -ijj. This is not too severe a restriction as the 
choice of function ip{ip) is always up to the designer. Performance characteristics of the system for this case are 
formulated in Proposition ^ 

Proposition 1 (Exponential Convergence) Let Assumptions\^\^hold and ip{ip) = Kip, K > 0. Then 
P6) function ip(x.{t),t) converges exponentially fast into the domain \ip{x{t),t)\ < O.5w|j0(O) — \\ff^£)p\-i- 



Furthermore, let Assumption^ hold, function a{x.,t) be locally bounded with respect to x and uniformly 



function q;(x, i) is persistently exciting: 

ft+L 

3L > 0, 5 > : j a(x(T), r)a(x(T), rfdr > SI \ft > 0, (15) 

where I G R'^^'^ - identity matrix. Then 

P7) both ip{x{t),t) and \\6 — \\ converge exponentially fast to the origin. 

It follows from Proposition ^ that if (p{ilj) = Kip then the estimate of the upper bound supj>f/ \ip{x{t),t\ as a 
function of t' for ip{'x.(t), t) in system H13|l exponentially converges into the domain determined by the parametric 

^Robustness and some other properties of algorithms 1121 are discussed in Appendix 1. 



Specifically, the following holds: \ip{:>i{t),t)\ < |V'(x(0), 0)|e--^* + 0.5./||^(0) - ^*||? 




bounded in t; for any bounded x there exi st Di > such that |z(x,0,t) - z(x,0 ,t)\ > Di\a{x,t)'^ {9 - 9 )\, 
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uncertainty, the values of controller parameters K, and adaptation gain F. Notice that this domain can be made 
arbitrary small, subject to the choice of the values of K and F. The rates of convergence are given by P6). In the 
case of persistent excitation an even stronger property is established. The system is shown to be exponentially 
stable with respect to the target manifold ■(/'(^i ^) = and point 6{t) = 6 . 

Despite properties of algorithms (|12|l such as improved transient performance of the closed loop system 
and their ability to deal with nonconvexly parameterized models, these algorithms are not realizable in the 
form of differential equations, as they depend on unknown parameters explicitly. It was proposed in |3U| 
to use special filters to estimate -0. While the approach of [30] is acceptable for systems with nonconvex 
parameterization, control system performance may be suboptimal due to estimation errors. The question is 
how to realize algorithms (|12|l in a form that depends neither on time-derivative "0 nor on its filtered estimate 
explicitly, nor on anything implying knowledge of unknown parameters 0. Our solution, as mentioned in Section 
1, is to use the finite form ||2Jl of adaptive algorithms instead of the differential form (|12|l . In the next section 
we study under what conditions algorithms (I12|) can be represented in finite form |(2J). 

3 Adaptive Algorithms in Finite Form 

The outline of the section is as follows. We start from a general case and formulate the conditions ensuring 
the realization of algorithm (|12|) in finite form explicitly, i.e., without any filters and further transformations 
of the closed-loop system. The conditions we impose involve the existence of the solutions of a system of 
partial differential equations. It is nontrivial to check these assumptions for nonlinear model ||2Jl. That they 
hold, however, can be demonstrated for some special combinations of plant models and goal functions i/'(x, f). 
Further, we consider extension of the proposed method to a broader class of nonlinear systems including systems 
with low-triangular structure. 

3.1 Explicit Realization 

Let us assume that in addition to the Assumptions ^{Sj that are sufficient for system (|13|l to have the properties 
Pl)-P7), the following hold 

Assumption 7 (Explicit realization condition) For the given functions a(x, i) and ip{x,t) there exists 
function 5'(x) such that the following hold: 



^{x) : a*(x,<)/ax2 = V'(x,t)(5a(x,t)/9x2) 



(16) 



Then realizations of the adaptive scheme described by equations (|12|) follow from the next theorem. 



Theorem 2 Let Assumption^hold. Then there is a finite-form realization of the algorithms ^1^) : 



e{x,t) 



F(6»p(x, t) + ei{t))- 0p(x, t) = V(x, i)a(x, t) - ^'(x, t) 



01 



.^(0(x, t))a{x, t) + a^'(x, t)/dt - V^x, t){da{x, t)/dt) 



(V'(x, i)Lf,a(x,t) - Lf,*(x, t)) - (0(x, t)Lg,a{x, t) - Lg,^'(x, t))w(x, 6>, t) 



(17) 
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Remark 1 It is easy to see from ()17|l and the theorem proof that reahzation of the algorithms 

= r(V; + /3(x,i))a(x,i), (18) 

where /3(x, t) is to guarantee at least the existence of solutions for the closed loop system, is also possible. 
Indeed, in order to realize these algorithms it is sufficient to replace equations for 6i in (|17|l by the following: 

Oi = /3(x,i)a(x,i) +9^'(x,t)/at-'0(x,t)(aa(x,i)/9t) - 

(^A(x,i)Lf,a(x,i) -Lf,*(x,t)) - (V(x,t)Xg,a(x,t) -ig,*(x,i))u(x,0,i), (19) 

One particular case of function (3{x,t) = (1 + (5(i))V'(x, t), 5 : i?+ 7?+, 5 G C" will be used later to show 
existence of the adaptive control algorithms for nonlinearly parameterized plants in the low-triangular form. 

Theorem |21 provides us with an answer to the question of existence of realizable algorithms that satisfy 
differential equations H12|) , thus ensuring the properties formulated in Theorem ^ ^^'^ Proposition ^ The 
disadvantage, however, is that the functions 'I'(x, i) in Assumption [7| are not easy to find. Existence of such 
functions itself is another nontrivial issue. For instance, if dimx2 = n and functions ijj{'x.,t), q:(x, t) do not 
depend explicitly on time t, then the necessary conditions for the function ^'(x) to exist is the symmetry of all 
matrices ^ ^■0(x) ^"g'x^'' ^ i i G {1, . . . ,n}. Nevertheless, despite difficulties in finding those functions ^'(x, t) 
that satisfy Assumption[71 there are several classes of dynamical systems with certain structural properties that 
immediately reduce Assumption 13 to more easily verifiable requirements. 

Corollary 1 (Single-dimension uncertainty-dependent partition) Let dim(x2) = 1 and function 
ijj{'x.,t)da{'x.,t)/dxn be Riemann-integrable with respect to a;„, i.e., the following integral exists 

*(x,0 = yV'(x,t)^^^dx„ (20) 

Then there is a finite-form realization of algorithms J^l°^) . 

Remark 2 Corollary allows us to turn the problem of searching for a function 5'(x, t) satisfying equation 
l)16|l into a problem of existence of the indefinite integral of a function with respect to a single scalar argument. 
It is clear from (|2(J|) that any one-dimensional system with integrable ^{x)^^^^ has a finite-form realization. 
An interesting example is the class of systems described by the following differential equations: 

Xi = /i(x) -f- .gi(x)u, i = l,...,n-l 

in = /„(x) +i/(x, 0) +g„(x)u, (21) 

where function /^(x, 6) satisfies Assumption|Sl which in turn is automatically satisfied if f(x, 6) linearly param- 
eterized or i^{x,6) = i/{x^6) and i'(-) is monotonic and belongs to a sector. In practice, the indefinite integral 
in H20() can also be replaced by ^'(x, t) — J^'^^q-^ ifii^, t)^^^^-^dxn. 

Equations of type (|21|l describe a class of dynamical systems in which the uncertainties are concentrated 
in a single equation. There are many mechanical systems described by equations H21() satisfying Assumption 
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{7\ (e.g., the simple classical equations for shaft dynamics with unknown load torque, rotating platforms and 
pendulums). For the case where the uncertainty is a single scalar and function a{x.) can be chosen as the goal 
function ^(x) = Q;(x), finite-form realization is also possible (^'(x) = ^a'^{x)). Another class of dynamical 
systems that automatically satisfy Assumption [7| is given by the following corollary. 

Corollary 2 (Independence on partition C2) Let function q;(x, i) be independent on C2, i- e., for any 

X2 G C2 the following holds: 9q;(x, i)/c'x2 — 9q:(xi ©X2,i)/9x2 — then there is a finite-form realization of 
algorithms ^T^. 

Corollary 121 conditions are equivalent to the fact that the plant dynamics can be described by system 

Xl = fl(xi ® X2) + gl(xi ® X2)W 

X2 = f2(xi®X2)+//(xi,0)+g2(xi©X2)u (22) 

and that dip{xi ©X2,i)/9x2 — A(xi,t), where A : _R™ x R ^ i?"^™ is a differentiable function with known 
derivative d\{xi,t)/dt. Therefore it is possible to derive from Corollary|21that every error model: tp = —(p{ip) + 
z(u){t),6) — z{uj{t),0), where Lj{t) : R — > i?", u; e is a function with known time-derivatives (jj{t), satisfies 
the sufficient conditions for realization of algorithm (I12|l in finite form. Indeed, this follows directly from 
Assumption [51 as functions a{x.,t) in this case are independent of x. Therefore if the derivatives a{t) are 
known, the finite form realization follows immediately from 

0(x,<) =^ r(0p(x,t) + 0/); 0p(x,i) = ^^(x,i)a(i) 
Oj = Lp{ip{x,t))a.{x,t) ~ ip{x,t)a.{t) 

This fact, along with decomposition (|22|l . will be used later in Section 3.2. 

So far, simplified conditions for the existence of the adaptive algorithms in finite form were derived from 
Theorem|21for those classes of nonlinear systems that have certain structural properties, such as single dimension 
uncertainty-dependent partition (Corollary ^ and equation (|21|) ') or independence of z(x, 0,<) on uncertainty- 
dependent partition X2 (Corollary [21). These structural properties allowed us to reduce Assumption [3 to inte- 
grability of a function with respect to a single scalar argument for a class of nonlinear systems. Taking these 
results into account, in the next section we present a technique that allows us to extend our method to a broader 
class of systems. 

3.2 Asymptotic Design via Embedding 

The main idea behind the extension of our results to a broader class of nonlinear systems is as follows. Instead of 
trying to find a general solution of equation H16() in Assumption[71 (which is a nontrivial task even if such solution 
exists), we transform the original equations into a form that satisfies much weaker requirements considered in 
Corollariesnand[21 This transformation should not necessarily be a one-to-one diffeomorphism, but the control 
goal reaching in the new state space should guarantee reaching the control goal of the original system. One 
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way to assure this is to embed the original system dynamics into one of a higher order, for which a finite form 
reahzation of the adaptive control algorithms is possible. 

Let us represent the partitioned system in the following way: 

Xl = fi(x) + gl(x)'U 

x^ = f^(x) + i.'(x,0)+g^(x)u; x'2' = f^'(x) + iy"(x,0)+g^'(x)u, (23) 

where X2 © X2 — X2, dimx2 = rrii, dimx2" = n — m — mi, < mi < n — m. Using the notations above, we 
introduce the following assumption 

Assumption 8 There exist 

1) a partition of the state vector x; x = Xi ® Xj © Xj , 

2) a system 0/ differential equations 

i = f5(x,ei); y? = he(0, (24) 
I e R\ : R"" X R'' X R+ ^ R\ e : R'' G C^; 

5j a function ^'(x, t) £ C""^, x = Xi © Xj © swc/i </iat the following conditions hold 

t)/dx'^ = ^(ic, i)(aa(i, i)/9x^) (25) 

Loo ^ ^ e Loo (26) 

/or any 6 E fig and t G -R+ along the solutions of the original system f^j. 

In addition to Assumption|Hl we would like to formulate two alternative assumptions which, if satisfied, will 
result in two different adaptation schemes with different performance and robustness properties. 

Assumption 9 Let system \24\l he given and 

z(x, e, t) - ^(i, e, t) G Loo] (V'(x, t) ~ V'(i, i))(aa(x, i)/9x'2)iv'(x, 6) G Loo- (27) 

along the solutions of (0), \24\l . 

Assumption 10 Let system \24\l be given and da.{x, t) / = 0, or tp{x,t) — tp[5(.,t). Furthermore, let 
z{'x.,6,t) — z{x.,6,t) G L2 for any E ilg, t > along the solutions of Q), i24}) - 

Sufficient conditions for the desired embedding follow from the next theorem. 

Theorem 3 (Embedding Theorem) Let function i/j{x,t) be given and Assumptions^]^\^ hold for system 
10). Then for the extended system 

X f (x) + i?(x, 0) + g(x)u 

i = fe(x,|,t); -h^(0, (28) 
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there exists control function m(x, hj, 6, t) 

u{x,h^,e,t) = (Lg(,)7A(x,t))-i(-^(^)-if^(x,i)-L^(-_e)V(i,i)-5V'(x,t)/ai) (29) 

and adaptation algorithms'^ : 0{i.,t) = T{6p{x,t) + Oi{t)), F > such that the following statements hold: 
P8) if \lp{i]-))\ > K\^\, K > and Assumption\^ holds then ■(/;(x, i), x, ^, £ Loo; 
P9) if Assumvtion \lUl holds then ip{'x., t) £ L2 Cl Loo, ^ £ -^2, ^(x, 6, t) — z(x, 6,t) £ L2, x, ^ G ioo/ 
if in addition derivatives 9'!/'(x, t)/9x, d'ip{:x.,t)/dt are uniformly bounded int and z{'x.,6,t) ~ z{Si,6,t) G L^o 

then -tjj e Loo, z{5<.,6,t) - z{x,6,t) e Loo, limj^oo V''(x(0, = 0. 

Theorem 13 states not only existence of the adaptive control algorithms but also provides us with exact 
equations for the adaptive control function. These equations are given by H29|l . 1)77(1 . which guarantee P8), and 
(|HTTl or ensuring P9) for da{5c)/dx'2 = or i/'(x) = V'(x) respectively. 

Although TheoremOguarantees reaching of the control goal and ascertains performance improvement (prop- 
erty P9)), it does not ensure the same properties of adaptive control algorithms as Theorem |21 does. On the 
other hand, the ability to deal with nonconvex parameterized systems is preserved, except for cases that do 
not satisfy Assumption |H1 The drawbacks of this narrower class of nonlinearly parameterized functions in the 
plant right-hand side and a slight degradation in performance are compensated by relaxing the requirement 
l)15|l of Assumption [7| Notice also that the difference in guaranteed performance reflected in P8) and P9) 
has the consequence that the dimensions of vectors are likely to be different in the both cases. Indeed, to 
ensure equality c)q;(x, t)/(9x2 = for arbitrary smooth function «(•), we must replace the whole vector X2 by 
hj(^). Therefore, in principle, embedding of the original system dynamics into one of a higher order is desired 
if improved performance and extended applicability are required. 

Theorem O offers a possible way to facilitate the search for function ^'(x, t) satisfying partial differential 
equation l|lt)|) as defined in Assumption [7| We replace the problem by one of searching for the embedding (|28|l 
which satisfies Assumption ISl and 1^ or [TUI The main obstacle, finding a solution to equation H16|) . is replaced 
with problem H25|l . the complexity^ of which should be reduced, as dimx2 < dimx2 if embedding into the 
higher-order dynamics is used. 

Indeed, according to Assumption |S1 and notations introduced above, the dynamics of the extended system 
can be described as 

Xl = fi(x) gi(x)w, X2' = f2'(x) + /^'(x, 0) + g2'(x)u 

= ^f5(x,|,i), i2" =f2"(x)+I^"(x,0)+g2"(xK (30) 

where vector Xi stands for the uncertainty-independent partition in the extended state space, and vector 
X2 is chosen to satisfy equation H25|) . Observe that function z{x,6, t) is independent of Xj and dimh^ = dimxj . 
Then for any h^: dimh^ > 0, we can conclude that dimXj < dimx2 = dimxj © X2. 

*The adaptation algorithms that guarantee properties P8) and P9) are given by equations I77i and 1811 . I82i respectively in 
Appendix 2. 

^Here reduced complexity means that the number of equations in the system is reduced. 
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Notice also that, by the appropriate choice of the dimensions of vectors ^ and (dimh^ = dimx2) in 
l)24|l . the dimension of vector can be reduced to unity. Alternatively, we may try to annihilate the partial 
derivative in (|25|) . Hence, eventually either Corollary ^ or Corollary [21 conditions will be satisfied for 

the extended system (|3U|I . This, in turn, implies that we can replace assumption (|25|l by a weaker requirement, 
such as integrability of the function with respect to a single scalar argument. 

After obtaining computable function ^'(x, i), the remaining problem is that we should be able to find 
an extension (|24l) that guarantees properties 1)2711 and (|26|l for the given partition x = xi © Xj © h^. If 
such an extension exists, then Assumption |S1 is automatically satisfied, and adaptive control algorithms follow 
immediately from Theorem |31 

Finding extension (|24ll that ensures boundcdness (and square integrability) of the differences z{x,6,t) — 
z(x, 0,t), (?/'(x, i) — V'(x, <))(9q;(x, t)/9x2)/>''(x, 0) is not an easy problem - taking into account that partition 
X2 is also uncertainty-dependent. It is possible to solve it using specially designed adaptive or high- gain auxiliary 
subsystems that track the reference signals X2 with the desired performance: z(x, 0,t) — z(x, 0, i), (?/'(x, <) — 
'0(x, i))^^|^^^i''(x, 0) G L2 n Loo- If, for example, partition X2 is linearly parameterized (i.e., v>"{x,9) = 
ri"(x)d), functions z{x,9,t), 'ip{x,t)^^^^r^i/'{:x.,6) arc locally Lipshitz in Xj and for any 9 £ Qg the follow- 
ing inequalities hold: |z(x,0,i) - z{i,e,t)\ < Ai(x, 4, i) ||x^' - h^(^)||, ||(V'(x,t) - V(x, ^)) ^"^^^/V ^x, 0) || < 
A2(x, ^, i) ||x2 — 115(^)11 then the suitable extension is defined by the following system: 

^1 = f^'(x)+r7"(x)^2 + A(x,ei)(x^'-^i)+g^'(x)« 

^2 = ri(x'2'-^i)V(x), ri>o, h^(0 = ^i, (31) 

where ^ = © ^2 ^^'^ M'^^ ^) — -^1(^1 ^) + -^2(^5 To show this, it is sufficient to consider the following 
Lyapunov's candidate: F(x,^) = O.5||(x^'-|^)||2-|-O.5||0-|2llr-i and observe that V < -A(x, |, i) ||x^' -^^ Ip < 
-(z(x,0,i) - z(i,6»,i))2 < 0. 

For a class of nonlinear systems with low-triangular structure 

Xn = fn{xi, . . . ,Xn,9„) + u + e{t), e{t) e L2, 9i eile (32) 

the suitable extension is guaranteed by Lemma 13 in Appendix 1. Then, combining the results formulated in 
Theorem 13 and Lemma |3 it is possible to show that our approach can be extended to a broad class of systems 
like those given by equations l|32(l . Let functions fi{-) in (|32|l satisfy the following assumption 

Assumption 11 Let there exist smooth functions Di{-) : W x x R ^ R such that for any 9i G the follow- 
ing holds: {fi{xi,. ..,Xi,9i)~ fi{x[, . . .,x'^,9i)f < 5f (x^, x^) [[x^ - x^p, x^ = (xi, . . . ,0;^)^, x^ = {x[, ...,x[). 

It is clear that Assumption 1111 holds for those functions /i(xi, . . . ,Xi,Oi) that are, for example, Lipshitz in x. 
The results for low-triangular systems H32|l are formulated in the next theorem 
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Theorem 4 (Finite Forms for Low- Triangular Systems) Let system if^jj)) and goal function ip{xi) = be 
given, and there exist functions ai{xi, . . . , Xi) such that Assumptions\^\^hold for the functions fi{xi, . . . ,Xi, 0i) 
in respectively. Furthermore, let fi{xi, . . . ,Xi, 0i) satisfy Assumvtion \l U a.i(xi, . . . , Xi), fi{xi, . . . ,Xi, 6i), 
i ^ 1, . . . ,n, -01 (xi) be smooth and the following condition holds: ip{xi) £ L^o ^ Xi £ Loo- 
Then there exist an auxiliary system 

k = f4(x,eiv), loGi?", /> = f,(x,|,iy), ivoGi?™, (33) 

as well as smooth functions ipi{xi, t), i — 1, . . . ,n, Op{'X., control u(x, 9, ^, u), and adaptation algorithm 

0(x,|,0j) = 7(^p(x,O + 0/), 7>0, bi^{gi^,e,^,iy), 

such that 

1) 'ipi{xi,t),Tl> e ^2 n Loo, ■ijj,Tpi £ L2, i = I, ■ ■ ■ ,n 

2) 6 £ Loo and u{-x, 0, u) - u(x, 0„, |, v) £ L2 

3) x,i,i^ £ Loo 

4) if e{t) £ Loo then ipi £ Loo, and limt^oo i'{xi{t)) = 0, linit^oo Aixii't),t) = 0, i 1, . . . , n. 

Theorem 0] extends the apphcabihty of algorithms in finite form to systems described by equation H32|) . 
Relying entirely on Lemma 13 (Appendix 1) and Theorem 13 Theorem 0] allows us to design adaptive control 
algorithms for cascades with nonlinear parameterization without the need for damping nonlinearities. However, 
performance is weaker. For instance, decrease (non-increase) of the term ||0 — 0(t) ||p_i is not guaranteed in this 
case. Nevertheless, adaptive control algorithms in finite form, in addition to their ability to deal with nonlinear 
parameterization, still guarantee certain improvements in performance. For instance, square integrability of the 
control effort due to adaptation (statement 2) of the theorem) and tpi{xi, t),'ipi £ L2 C] Loo are ensured. In the 
next section we illustrate our method with the examples. 

4 Examples 

Let us consider the following system: 

xi ^ xl8o + x2; X2 ^ xi9i + X2O2 + u, (34) 

where parameters 6o,0i and 02 are assumed to be unknown. The control goal is to steer the system towards 
the following manifold: a;i — 1 = 0. To design adaptive algorithms in finite form for system (|34|l . we follow the 
steps of Theorem 0] proof: 

1) Intermediate control design. Derive control function ui(a;i,6'o) such that for the reduced system 

xi^xleo + ui[xi,9a) + ei{t), ei{t) £ L2; = Oo^xi) + O^jit) 

reaching of the control goal is guaranteed: '4){xi{t)) — xi{t) — \ Q as t ^ 00. Moreover, function 
ui(xi,6q{xi,6qj)) should ensure that ip,ip £ L2. 
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2) Embedding. Extend the system dynamics (or embed it into) with auxihary system 

i = fd^,^.v)] V = f^{-^,^,v) (35) 

in order to guarantee that 

M(a:i,^o(a:i,^o,/)) -w(^,^o(C,^o,/)) e L2, xi - ^ e L2 (36) 

3) Control function design. Introduce new goal function '02(a;2, t) = ^2 — iti(^, 0o(Cj ^0,1)) a-nd derive control 
function u{xi,X2,^,t) such that -02 G L2, ip2 ^ L2 f] Loo- The last automatically implies that xi — x^Oq + 
X2 = x\6o + ui{xi,eii{xi,9o.i)) + n{t), where ^(t) ^ X2 - ui{xi,0a{xi,§oj)) = {x2 ~ ui(C, ^o(^, ^0,/))) + 
(wi(^, 6*0(^1 ^0,/)) — ui(xi, 6*0(2:^1, ^'o,/))) G -^2- Therefore, according to the choice of function ui{xi, Oq{xi, ^o,/)), 
control u{xi,X2,S.,t) guarantees that ^j(xi{t)) — > as t — > 00, tp,ip E L2. 

We begin by determining the function u{xi, Oq{xi, Oqj)). Let ui{xi, Oq) = —{xi — l) — 9oXi, where satisfies 
the following differential equation: 

00 = {xi - 1 + xi)xl (37) 

It follows from Lemma |21 that control function ui{xi,9o) with algorithm (|37|l guarantee that Tp^ip G L2, 
tp{xi{t)) as t — > 00. According to Theorem |3 finite form realization of H37() can be given as follows: 
6'o(a;i, 6'o,7(t)) = l/Sxf + 6oj{t); Oqj = {xi — l)x\. Substituting this into ui{xi,9q) we get the following 
expression for ui{-): 

ui[xi,0o{xi,6o.i)) = -(.Ti - 1) - 1/3x1 - x\9oj[ty, §oj = -tlj{xi)ai{xi) = {xi - l)xl. (38) 
Thus step 1 is completed. 

Let us design system H35|l which guarantees that 136|l holds for function H38(l . First consider the difference: 

uixi,9oixi,9o.i)) - uit 9oi^, 9o.i)) - -{xi - 0(1 + [xi + Oki + + 4^ + + ^if + (39) 

and denote i^(a;i, ^,^7, 0) = {l + {xi+£,)9oj + l{xf + xl(,+xlS,'^ + xi£,^+£,'^)). It follows from Lemma|2|that there 
exists system (|35|) such that condition H3t)|l holds. In fact, this system can be given by the following equation 

C = {xi-^)iF^ixi,tOo,i) + ^)+49^ + X2 (40) 

where 9^ satisfies the following differential equation 9^ — [xi — ^ + ii — £,)4- Fhiitc form realization of this 
algorithm^ follows from Theorem [3 and it can be written as: 

9^ = 1/3x1 + 9^j; ki = i^i-04~4ii^i'-OiF^i^uCOoj) + l)+49^ + X2) (41) 

"introduction of algorithms I41i is not necessary here because the original system is linearly parameterized, and condition 1361 
can be satisfied even with conventional (gradient) adaptation schemes. Nevertheless, we would like to keep the consistency of our 
current calculations with those steps made in the proof of Theorem |1] in order to illustrate what would happen if the right hand 
sides are nonlinearly parameterized. 
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Taking into account (|41|l and (|4U|) system H35f) which ensures can be represented as follows 

ki = ixi~Oxl~xl{{xi-OiFHxi,^,0oj) + l) + lxl + 9^j{t)xl+X2). (42) 

Therefore, step 2 is completed as well. To conclude the controller design let us consider new target manifold 
X2 - ui{^,0o{^,9oj)) ^ and goal function 'ip2{x2,t) = X2 - ui{£,,6o{(,,9oj)) = a;2 + C - 1 + + ^o,/C^- Let 
us write function '02(') derivative with respect to time t: 

= X2 J- ^ -T 0QJ ^Xiei+X202 +U + ^ {xi-l)X^ + 

800.1 

(1 + 1^ + 2^0oM^i - OiF'ixutdoj) + 1) + + 0u{t)xl + X2) 
Therefore, control function 

U = -i0l-X202 -^''{xi-l)xl^{x2+^-l + ^S,'' +0O,lf)~ 

(1 + + 2i0oMxi - i){F\x^,^, 0o,i) + 1) + ^1 + 0u{t)xl + X2) (43) 

results in the following error model: "02 = —'>p2{x2, t) + xi0i + X202 ~ xi0i — a;202- Taking into account condition 
we can rewrite derivative ?/'2 as -02 = —^p2{x2,t)+^0l+X202—^0l—X202+£{t), where £(i) ~ ixi—^)0i G L2- 
It follows from Lemma |21 that adaptation algorithm 

^1 (V'2(a;2,i) +'02)q!i(C); ^2 = (02(2:2,0 + ^2)0^2 (2:2), ai(0=^, Q;2(a;2) = a:2 
guarantees 02 G -^2 H Loo and -02 G ^2- Realization of algorithms 144|) can be obtained from Theorem |21 

ei{x2,(Joj,t) = {x2 + i-i + ^e + 0oA^K + ki(ty^ ki = ix2 + ^-i + le + Oo,i^^K^-0 

02ix2,Cjo.I,t) = -^+02j{t)- 02j^{x2+^-l + -e +kl(')x2 + ^i+-^eo.I. (44) 

where ^2{x2,tOoj) = J Mx2,t)^^^dx2 = f + (^ - 1 + + ^o,,e')a;2. 

We would also like to compare performance of the proposed adaptation scheme with adaptive backstepping 
control algorithms. Adaptive backstepping design for system H34|l according to 8 results in control algorithm: 

Ui = —2X2 — {Xl — 1) — 03Xi — xf{xi — 1) — 20^X1X2 — (xf + 20^x1)0 — Xi0i — X202 

= (x2+a;i-l + 42;?)a;i(l + 2^3x1); ^1 = (2:2 + xi - 1 + ^3a;i)a;i 

§2 = (X2 + Xl - 1 + ^3x^)2:2; ^3 = {xi - l)xl (45) 
Adaptive backstepping with tuning functions results in 

Ui = -{X2+Xi-1+Xj0)-{XI-1)-{1 + 2X10){X2+0XI)-XIT-Xi0i-X202 

^ t; t ^ {xi ~ l)xl + {x2 + Xl - 1 + xl0)xl{l + 2x10) 

01 = (a;2 +a;i - 1 + a;?^)xi; ^2 = (a;2 + xi - 1 + x?^)a;2 (46) 
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Figure 1: Plots of system l|T7|l trajectories with control functions (thick sohd lines), (dotted 

line), l|46(l (dash-dotted line), and (|48|l (thin solid line). Plot a - xi as a function of time, 6 - 0:2 as a function 
of time, c - as a function of time, c? - u as a function of time. 



We simulated the adaptive system dynamics for the following set of parameters and initial conditions a::i(0) = 
2, 0:2(0) = 0.2, 61 = 1, 613(0) 61(0) = 3, 6li(0) = 6(2(0) = -2, ^(0) = 0, 6(0) = 0, k = 10. Initial conditions 
for ^1,7(0), 02,/(O) and ^3,7(0) where chosen to satisfy ^i(O) = ^2(0) = -2, ^3(0) = 3. Parameters 6*1, 6*2: 0i = 1, 
62 = 0.5. As an additional measure of performance, we introduced the variable A9 which indicates the distance 
in the controller parameter space between the estimates and real values of the parameters. Simulation results 
are presented in Figure 1. In Figure 1 thick solid lines show the system dynamics with algorithm (|43() . H44() . 
dotted lines show the system behavior with algorithm (|45|l . and dash-dotted lines correspond to algorithm H46|) . 
It turns out that system H34(l with algorithm 1)46(1 also reaches the goal manifold, but after around 400 seconds of 
modeling time. We can see again that transient performance of the adaptive algorithms proposed in the paper 
is better than that of conventional algorithms. In addition, we calculated the integral / — u\{T)dT, T — 500, 
for every controller along the system solutions. The values of the functional / indicate how much energy is spent 
to achieve the control goal. For control function (|43|l . (|44() / = 627.10, for adaptive backstepping controller (|45|l 
/ 13329.28, for controller ((46|l / = 263872.58. To illustrate the ability of our algorithms to deal with nonlinear 
parameterization, we change H34|) to 

±1 — x\9 + X2] ±2 = 5tanh(a;i6'i + .X2^?2) + "2 (47) 
Nonlinearity tanh(a;i0i+X202) satisfies Assumption|31with respect to function q;(x) = (xi, x^)^ and, in addition. 
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Assumption El is also satisfied for any bounded xi and X2- Then according to Theorem^ control function 

(1 + + 2^0oM^i ~ i){F\x^,i, §0,1) + 1) + + e^,i{t)xl + X2) (48) 

along with 14411 guarantees that V'l ! '02 , '01 , V'2 S L2 H Loo ■ The simulation results of system (|47ll with control 
algorithm (|48|) . H44|) are given in Figure 1 (thin solid lines). The value of functional / for this case is 3186.83. 

5 Conclusions 

The method proposed in this paper suggests a new methodology to design adaptive control algorithms. Our 
method to design of the adaptation schemes is consistent with recent trends in adaptive control, for instance, [12j, 
where nonlinear controllers are proposed to adaptively stabilize linear plants. Indeed, when derived for linear 
systems algorithms in finite form will also result in nonlinearities in the controller. These nonlinearities are to be 
introduced, in particular, to improve the performance of the adaptive system. In contrast to jl2| . we show not 
only that the L2 and L^c norm bounds are computable for the state vector, but also that properties Pl)-P7) are 
ensured. The method, however, is different from conventional approaches, as it is not restricted by realizability 
issues. While in conventional parametric adaptive control the realizability of adaptation schemes in differential 
form determines the properties of the resulting systems (including poor performance and restricted applicability), 
in our method we first determine the desired properties of the controller (Theorem^ Proposition^and Lemmas 
^-l^l and only then deal with the realizability problem. In order to realize the adaptive algorithms in finite 
form explicitly, i.e. without extension of the system state space, special restrictions formulated in Assumption 
[7|are to be satisfied (Theorem |21 and Corollaries ^ IS ■ 

To realize the adaptation algorithms that do not satisfy the explicit realizability conditions formulated in 
Assumption we embed the original system into a system of higher order. This system should satisfy a-priori 
certain structural conditions that are formulated in Corollariesnand|21 These two ideas (design of an algorithm 
aiming for its best properties, not its realizability, and design of an embedding for realization) result in a new 
method, which is shown to be applicable to a sufRciently large class of systems with nonlinear parameterization, 
e.g., like those given by systems (|32|l . It is very important that no damping or discontinuities are injected 
directly into the control function in contrast to [3 115L . 

In the present article we hope to have extended the scope of applicability and performance of adaptive 
control algorithms. Our results to date are applicable to the full-state feedback case only. Extension of the 
results to the output-feedback case remains a future study topic. 
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6 Appendix 1 

In this section we consider auxiliary and technical results that are used in the paper. Let the system dynamics 
with respect to the function ip{x,t) be described as follows: 

= -ip{tp) -\- z{yi, e, t) - z{x, e, t) + e{t), (49) 
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where function e : R+ R, s E C'^ models unknown disturbances due to unmodeled dynamics or measurement 
errors. In addition, we will assume that the adaptation algorithms are affected by a disturbance: 



Lemma 1 Let error model be given, S,e G Lao, W{'>t-')\ > K > and Assumptions^\^^ hold for 

e = 0. Then ?/'(x, i) and 6, x(t) are bounded for the error model \4-y\j with algorithm 

e = (r(?/> + ^(V'(x,i)))a(x,i) +(5(t) - A > 0. (50) 

Lemma^proof. Denote Q{^p) — ip{<^)d(^ and consider the following function 

v{i:,e,e*) = 2{D - Di)Q{^) + o.5\\e - e*\\l-^. (51) 

Its derivative satisfies the following (z{x, 9, t) — z(x, 9, t) + e{t) = (p{ijj) + ijj due to equation (|49|l 'l: 

V = -2{D-Di)ip{il;)i^ + {e-9*f{{ip{il;)+i;)a{x,t)+S{t) + X9) = -2{D-Di)ifi;- 
(z(x, e, t) + z(x, e, t)){e - e*fa{x, t) + e{t){9 - 9*foL{K, t) + {9- e*f{S{t) - A^) 

From Assumptions [S] El it follows that 

-(z(x, 9, t) - z(x, 9, t)){9 - 9*fa{ii, t) + e{t){9 - 9*fa{x, t) - £l^l!£t + Rl^l^ < 

< -D{zi^, 9, t) - z(x, 9, t) f + i?i|£(t)||z(x, 9, t) - z(x, 9, t)\ - ^^^^ + ^^^^ 

= -(i?-i?i)(z(x,0,i)-z(x,0,i))2-i?i (^|z(x,0,t)-z(x,0,t)|-^^ ^52) 

Then 

mi 



V < 2{D - Di)if{iP)iP - [D - Di)(z{-K,9,t) - z{-K,9,t)Y - Di |z(x, 6», i) - z(x, 6/, i)| - 



2 



-^lli -{9 9 f{S{t) + X9) < -{D - D,){^^{^) + - 2{^{i^) + i^)e{t) + e^{t)) + 



4 



DisHt) 

— A\\0 — 



< -{D -D,)[l-^\ W + + 2{D - D,)AleHt) + - A fl - ) ||0 - ^ 



^2 ^i;^ . (^■i) 



where Ai, A2 > 1. In the lemma conditions (5, e G Loo- Therefore, taking into account estimate (|53|l and 
inequality |</'('0)l > ^IV'L we conclude that derivative V is negative-definite for any tp,9 that belong to the 
following set: 



nt>o^U,9 



4 / ^ 4A 
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Hence 6, ijj{x.,t) are bounded. T/ie lemma is proven. 

Let £ e L2 and 6{t) = 0. In this case it is possible to show that the control goal is reached in the closed loop 
system with slightly modified version of algorithm H12|l . 

Lemma 2 Let the following error model be given 

^ = -ifWa + F{t)) + z(x, e, t) - z(x, e, t) + £(t), (54) 

where F : R+ ^ R+, F{t) £ C°,, e(t) e C°, e(t) G L2, Assumptions \MB hold for F{t) = 0, e{t) = and 
adaptation algorithm satisfy equation 

b^T{lP + <f{^P){^+F{t))cx{x,t). (55) 

Then 

1) V'(X, t) e Loo, ¥'(V'(X, t)) e i2 n Loo, v/^F(^V5(V'(x, t)) e L2; e Loo 

2) z{x,e,t) - z{x,e,t) e L2. 

If F{t) e Loo then 

3) i^€ L2. 

If in addition functions e(t) £ Loo and function z{-x.,6,t) is locally hounded with respect to x, 0, uniformly 
hounded with respect to t then 

4) ■(/'(x, t) ^ as t —^ 00. 

Lemma\^proof. Function e{t) G L2, therefore integral e^{T)dT < 00. Consider the following function 

T/^(0,r,O = ^ e\T)dT + -\\9- 6*111^,. 

Its time-derivative can be written as follows: 

- -l/4D,e\t) + iip{i;)il+F{t)) + ^)ie-0fa{^,t) 

= -l/4:Die^{t) + (z(x, e, t) - z(x, 0, t) + e{t)){e - 0*)^a(x, t) 

Taking into account inequality H52|) we can write the following estimate for Vf. 

< -{D - D,){z{^, e, t) - z(x, 6, t)f = -{D - Di){^{^:){l + F{t)) + ^ - e{i)f (56) 

It follows from that zfx. 0. t) - z(x. e*.t) e L2. Let us denote ^(t) ^ e{t) + z{x{t),e{t),t)- z{x{t),e* ,t). 
Taking this equality into account, error model 154|) can be written as follows ip — — (1 + F{t))(p{ip) + fi{t), where 
function fi{t) G L2 as a sum of the functions from L2- Consider the following nonnegative function Vi{'ip,t): 

V,{^,t)= + f^\r)dT 

Its time-derivative is: 

V, = -vmi + F{t))v{ij) + ^(^)Ai(t) - \ii\t) ~ Fit)ip^{ij) - (^^(V) - i^w) (57) 

22 



It follows from inequality (|57|) that -0(x, t), (^(-0(x(t), i)) £ L^o- Furthermore, F(t)Lp{'4>{'x.{t),t)) G L2 and 
((y5('i/'(x(i), t)) — fi{t)/2) e L2- Given that /^(t) G L2 it is clear that (p{tp{x{t),t)) G ^2- Hence statements 1) 
and 2) of the lemma are proven. Let F{t) G £00 then (1 + F{t))Lp{ip{x{t),t)) G £2 and therefore according to 
H54|l ^ £ L2 'AS well. Thus statement 3) is proven. To show that 4) holds it is sufficient to notice that is 
bounded due to H56() . According to Assumption Estate x is bounded as t) is bounded. Then -0 is bounded 
if e{t) is bounded and function z(x, 9, t) is locally bounded. Hence applying Barbalat's lemma we conclude that 
^/'(x, t) ^ as i 00. The lemma is proven. 

Lemma 3 Let system: 

ii = fi{xi, . . . ,Xi,9i) + f3i{-x-,t), i = l,...,n (58) 

and smooth function u{x,z,9q) : i?" x i?™ x R"^ R, be given. Let us assume that 9q G Qq, Qq be bounded 
and there exist smooth functions i^(x, x', z), Di{x., x'), i — 1, . . . ,n such that the following properties hold: 

(u(x,z,6>o) -u(x',z,6»o))^ < \\x-^'\\^F^{ji,x',z) V^o e l^o,x,x' G i?" 

{f,{xi,...,x,,9,)- f,{x[,...,xl9,)f < ||i,-x^f^2(x,,x^ Gr!e,x„x^ Gi?" 

Xi = {xi,...,Xi,0,...,0)'^, X- = (x'l, . . . ,a;-,0, . . . ,0)^ 

Lei us also assume that there exist Q:i(x) such that Assumptions\^\^hold for the functions fi{xi, . . . ,Xi, 9i) 
respectively. 

Then there exist ^(t) : R ^ R" , i'{t) : R ^ i?™, smooth functions f^(-)j fj^(') o,nd corresponding system: 

«> = f,(x,|,z,iv), z^oei?'" (59) 

SMc/i that 

1) u(x, X, 0o) " u(qi, z, ^0) e £2, i = 1, . . . , n 

/i(a;i, . . .,Xi,9i) - ftiCi,. . . ,£,i-i,Xi,9i) G £2, i 2, . . . ,n; 

5; X G £00 ^ ^,1^ e £00- 

Lemma\^ proof . For the sake of notational convenience we would like to use the following notations: 

fi{xi, . . . ,Xi,9i) = fi{x.,9i), V'Ci ^Xi-£,i 

£i{t) fi{qi-2,di) ~ fi{q.i-i,di), i^2,...,n. (60) 
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Consider the following system of differential equations: 

k 

6 = ((F/(q,_i,q„z)+^5|+i(q,_i,qO) + l)(a:,-C»)+/.(q.:-i,%)+ft(x,i) (61) 

k 

% = 7?.('A?.((-Fi^(q»-i,qj,z) +^l)|+i(qj_i,qj)) + 1) + ^/'^Jaj(q»-i), 7?. > 0, i = 
where Fi(q,;_i, q,, z) = F(qi_i, q^, z). Taking into account H58|l and let us write the following error model: 

k 

k 

It is clear that trajectories of system for fc = 1 satisfy the following condition u(x, z, 0o) ^ '"(qij z, 9q) G L2 
(this follows directly from Lemma [21). Consider the case when k — 2. Taking into account the equations for 
and 6^^ we can derive from Lemma |5| that 

(qo,qi,z) + 5|(qo,qi)(a;i - ^1) G L2, - 6 G -^00, £ ^oo 

Hence ^(qo, qi, z)(a;i -^1) G L2, L'2(qo, qi)(a;i - Ci) e L2 as (qo , qi , z) + L»|(qo , qi ) > |Fi(qo,qi,z) 
and \/ (qo , qi , z) + -D| (qo , qi ) > |-D2(qo, qi)|- Therefore, we can conclude that 

u(qo, z, Oo) - M(qi, z, Oq) G L2, /2(qo, ^2) - /2(qi, ^2) e L2 

as |M(qo,z,0o)-'"(qi,z,0o)| < |-F'i(qo,qi,z)(a;i-6)l e ^2, |/2(qo, 02)-/2(qi, ^2)! < |52(qo, qi)(2;i-6)l e -^2- 

Notice that £2(0 = /2(qo, ^2) — /2(qi, ^2) = /2(x, ^2) — /2(qi, 02)- Therefore, the equations for ^p^^ become as 
follows: 

= -(F|(qi,q2,z) + l)V^, -/2(qi,%)+/2(qi,02)+£2(t) (63) 
% = 76(V'C2(A^(qi'q2:z) + 1) + V^?2)"2(qi), 742 > 0> e(^) ^ ^2 

Hence, applying Lemma El to system and taking into account that w(qi,z, 0o) — '"(q2,z,0o)| < 

|-F'2(qi, q2, z)(a;2 — ^2)! we can conclude that 

u(qo,z, 0o) - u{qi,z,9o) G L2, a;i - G Loo, e ioo 
u(qi,z, 0o) - w(q2,z, 0o) e -^2, a;2 - 6 e Loo, % e Loo 

/2(qo, ^2) - /2(qi, ^2) e L2, 

and, subsequently, u(qo, z, 69) ^ w(q2, z, 6q) G L2 as a sum of two signals from ^2- 

Let us now consider arbitrary 2 < k < n. It follows from Lemma (21 that for the error model with respect to 
function ipi^-^ : 

k 

= ~((Ff(qo,qi,z)+^^-+i(qo,qi)) + l)^6 -/i(qo,%) + /i(x,ei) 



24 



and corresponding subsystem 

k 

hi = 76(V^«i((A'(qo,qi,z) + 5I^I+i(qo,qi)) + 1) + V'^Jai(qo), 



one can derive that y Fj^(qo, qi, z) + J2j=i ^|+i(qo, qi)(a^i — ^i) G -^2- This in consequence imphes that 

m(x, z, Oq) - u(qi, z, Oq) G L2 
/,(x,6>,) - /,(qi,0,;) eia ,i = 2,...,k 

Hence we can write the error model for t^^^ in l(j2|l in the foUowing form 

k 

V-^ = - ( (^1 (qi , q2 , z) + ^ 5'+ 1 (qi , q2 ) ) + 1 - /2 (qi , ) + /2 (qi , 02 ) + £2 (<) 

i=3 

where £2{t) G i2- It foUows from Lemma |2 that 

ii(qi, z, Go) - ii(q2, z, 0o) e -^2 

/i(qi,0O - /i(q2,0z) e L2 ,i = 3, . . . , /c 

in system (|61|) . Notice also that /3(x, 9^) — /3(q2, ^3) G L2 as a sum of two functions from L2- By the similar 
reasoning it is can be shown that for any 2 < i < n we can represent the error model system (|62|l as follows 

k 

4 = -((-F;'(q.-i,q.>z)+ ^ ^f+i(q,_i,q,0)) + i)V«, -/^(q»-i,^«J + /^(q»-i,0O + £.W, 

j=i+i 

where £i(t) G L2. Therefore, using Lemma |5| again we can conclude that 

M(qj_i,z,0o) - '"(qj ,z, ^o) G ^2 
/i(qj-i,0i) - /j(qi>0i) G ^2, aJi-^iGioo, % G Loo, 

The last, however, implies that m(x, z, 0o) ~ ""(qi, z, 6q) G L2- In order to complete the proof we have to make 
sure that system (|61f) is physically realizable. In particular, realization of subsystems 

k 

% =7«,(V'6((^'(q»-i,q»,z)+X!^l+i(^»-i'^»)) + l)+4')"»(^*-i)' 74, >0, 1 = 1,. ..,fc, (64) 

shell not be dependent on any uncertainties Oi. It follows, however, from Theorem |5] that there are realizations 
of algorithms (|64|l in finite form: 

%(q,_i,6,t) = 7^^(%,p(q,_i,e,) + %j(t)), >0 

^«.,p(q»-i,6) = %(a;j,6)aj(qi-i) - *«,(qj-i,Cj) 

= ((^^(q^-l,q^,z) + + + 

^a%(q,_i,6)^ , , > ^9Q;i(qi_l) : 

L - L ' ) 
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Notice also that if x G Loo then ^ e Loo and hence 6^i,p{(li-i,£,i) G Loo as 0^i.p{cii~i,£,i) is smooth. Given 
that 0^. — j^.{0^.^p{qi-i,£_i) +0^.j) and both 0^.,0^.^p{qi-i,S^i) e Loo then we can conclude that O^.j e Loo 
for X G Loo. 

It is easy to see that denoting 1/ = ^^.j we can transform system (|61|) . (|65f) into (|59|l which satisfies 
statements l)-3) of the lemma. The lemma is proven. 

7 Appendix 2 

Theorem^ proof. Let us consider the following positive-definite function: Vg{d,9 ) = ^\\6 — 9 ||p_i. Its 
time-derivative according to equations (|13l) can be derived as follows: ¥^{9,6 ) = {ip{il;) + ij){9 - 9 )^a(x,<). 
According to Assumption and equahty © it is easy to see that 

Vg{9,9*) = ~{z(x,9,t) ~ z{x,9,t)){9 ~ 9fa.{x,t) < ~D{z{x,9,t) - z(x,9,t))'^ ^ -D{if{ij)+^f < (66) 

Therefore Vg is non-increasing (property P2) is proven). Furthermore, integration of Vg with respect to time 
results in 

V^{9{0),9*) - V^{9(t),9*) > D fwr) + ^{^P{T))fdT > 0. 

Jo 

Function Vg is non-increasing and bounded from below as Vg > 0, therefore 

D f\i;{T) + ^(Mm'dT < V^{m,9) < oo. 
Jo 

Hence ((^(V') + ip) = (^(x, 9, t) - z(x, 0, t)) = (z(x, 9* ,t) - z(x, 0, t)) G La (property P3)). 

To prove property PI) let us consider the following fimction: V{ifi,9,9 ) = 2DQ{%1)) + V^{9,9 ), where 
Q^tp) = Jq ip{'^)d(,. Function V{ip, 9) is positive-definite with respect to '0(x, t) and 9 9 because of Assumption 
21 Its time-derivative obeys inequality: V{^, 9, 9*) < 2DLp{^)"^ - D{'iij + (pii'))'^ = -Dtp^{^) - Dip^ < 0. 

Therefore, function Viip, 9,9 ) is bounded and non-increasing. Furthermore 

oo>y(V'(x(o),o),0(o),0*)>\^(V^(x(o),o),0(o),r)-y(vXx(<),i),0(i),r) > of ^\ij{^{T),T))dT >o 

Jo 

cx)>l/(V'(x(O),O),0(O),0*) >y(V'(x(O),O),0(O),0*)-y(V'(x(t),i),0(i),0*) > D [ ij^{T)dT>0. (67) 

Jo 

or, equivalently, S -^2, v('0(^)) G ^2- Hence, property PI) is proven as well. The L2 norm bounds p4|l for 
ip{ijj) and -ip follow immediately from inequality (|67|l : 

||<^(V')ll2 < i?- V(^(x(o), 0), 0(0), 0*), ||^||2 < i?-iy(^(x(o), 0), 0(0), r ) 

The Loo norm bound for V(x(t),i) results from the inequality: V{ip{ii{0),0),9{0),9*) -V{ip{Ti{t),t),9{t),9*) > 
0. Consider function A defined as A(d) = maxi^idi/)! | J^^^^ ip{<;)d(^ — d} and notice that it is monotonic 
and nondecreasing. Therefore, given that /q (y9(i;)(i^ < j-^V{tp{x{0),0),9{0),9 ) we can conclude that 

IV'I < A ('2^y('0(x(O),O),0(O),0*)). To prove property P4) notice that function V{ilj{x{t),t),9{t),9*) is 
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bounded. Hence by Assumption 0] function ■tp{x{t),t) is bounded as well. According to Assumption ^bounded- 
ness of V'(x(t), t) implies boundedness of the sate x. In addition it is assumed that z(x, 6, t) is locally bounded 
with respect to x, and uniformly bounded in t. Therefore the difference z{'x.,6,t) — z{'x.,6,t) is bounded. 
Furthermore, by Assumption ^function (p{ip) G and therefore it is bounded as well given that ip is bounded. 
Hence ip is bounded and by applying Barbalat's lemma one can show that il){n.{t)^t) — > at t ^ oo. 

To compete the proof of the theorem consider the difference z(x, 0,t) — z{x,6,t). Let function (p G C^, 
function z(x, is differentiable in x, 9; derivative dz{x,6,t)/dt is bounded uniformly in t: function q;(x, i) 
is locally bounded with respect to x and uniformly boimded with respect to t, then d/dt{z{'x, 6, t) — z(x, 9, t)) 
is bounded. On the over hand there exists the following limit 

lim / (z(x,0,t)-z(x,0,t))2= / (z(x,0,r)-z(x,0,r))2 < -F^(0(O),0*) 

as Jg (z(x, 9, t) — z(x, 9, r))^ is non-decreasing and bounded from above. Hence by Barbalat's lemma it follows 
that z(x, 9, t) — z(x, ^, r) ^ as i ^ oo. Notice also that ip{x{t),t) — » as t — > oo. Then ip ^ a.s t ^ oo. 
The theorem is proven. 

Proof of Proposition^ Consider the following integral'' /q (V'(t) + tp{'il:{T))'^ dr . It was shown in Theorem 
n proof that /q('0(t) + (p{'ijj{T))^dT < jjj\\9{0) — 9 ||p_i along system ((T^ solutions. Let us define fi{t) = 
ip{t) + (p{ip{t)). In the other words 

= (68) 

where ii'^{T)dT < jjj\\9{0) — 9 ||p_i. According to the proposition conditions, (p{i^) — Kip, it is possible to 
derive the solution of equation as follows %p{t) — ■(/'(0)e~^* + e^-^'-*^'^^ ^{T)dT. Hence 



IV'Wl < \^iO)\e-'^' + \l { t e-^(*-^)ti{T)dT] < |V(0)|e-^* + W /*e-2^(*--)dr /V'W^^^ 



/ V •'0 Jo 



< \m\e-''' + y^\\m~e*\\l_,. (69) 
Property P6) is thus proven. In order to prove property P7) consider 

9^T{^ + V5(V'))a(x, t) = r(z(x, 9, t) - z(x, 6>, t))a(x, t). 

Function 

Di\a.{yi,t)'^ [9 - 9*)\ < \z{:>i,9,t) - z{x,9,t))\ < D\cx{x, t)'^ {9 - 9*)\ 
a{y:,tf{9-9*){z{yi,9,t) - z{y:,9*,t)) > V z{x,9*,t) ^ z{x,9,t). 

Therefore, there exists Di < K{t) < D such that 

9 = -K{t)Ta{x,t)'^{9 - 9 ) a.{x,t) = —K{t)ra.{x,t)a.{x, tf{9~9). 

'^That we substitute the arguments of the functions i/>(-) and with t means that we consider them as functions of time. 
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Hence 

9{t) -9* ^ e'^'I'o ''(^)«(x(r),r)«(x(r),r)^dr^^^Q) _ g*^ ^jq^ 

Consider the integral T K{T)a.{x{T),T)a.{x{T),T)'^ dr for t>L 

rf K{T)a{-x{T),T)a{x{T),T)'^dT >TDi f a{-K{T),T)a{x{T),T)'^dT, 
Jo Jo 

where a('x.(t),t) is persistently exciting. For any t > L there exists integer n > such that t — nL + r, 
r e R,0 < r < L. Therefore 



TDi I cx{yi{T),T)cx{x{T),TfdT>TDinSI> 
Jo 

Then taking into account ()70|l one can write 



t-I 



||^(t)-^*ll<lle(-^*+^)||||0(O)-r||, (71) 

i. e. 6{t) converges to 6 exponentially fast. It means that there exist positive constants \ > 0, \ ^ K and 
Dg > such that \\e{t) - < e^-^*||0(O) - 6*\\Dg. It follows from Theorem [T] that ?A(x(i),i) is bounded. 
In addition due to Assumption ^ we can conclude that x is bounded as well. By the proposition assumptions 
function q;(x, t) is locally bounded with respect to x and uniformly bounded in t. Therefore, there exists Da > 
such that \a{x,t)'^{e{t) - 9*)] < Da\\9{t) - 9*\\. Taking into account that z{x,9,t) - z{x,9,t) = and 
\z{x, 9, t) — z(x, 9,t)\ < D\a{x, t)'^{9{t) — 9 )\ one can derive from (|68|l the following estimate 

\m\ < + \\m - 9*110^0^0 e-^'^'-^h-'^dT < |V(0)|e-^* + ^^^^^ \\m - (72) 

The proposition is proven. 

Proof of Theorem\^ The theorem proof is quite straightforward and follows from explicit differentiation of 
function 0(x, t) with respect to time: 9{x, t) = T{9p + 9i) = T{i>a.{yi, t) + ■(/'q;(x, t) - \i'(x, t) + 9i). Notice that 

V-a x, t) - * x, t + ^/ = X, i) ^ ^ xi + 1/; x, t)^^^2 + X, t) \: ' - 

0x1 0x2 at 

g^(x,t) . d^{x,t) . d^{x,t) , 

— 5 xi X2 1- 9i (73) 

axi 0x2 ot 

According to Assumption [7| ^^q^^*'^ — ^(x, t) ^"'g^f' ■ Then taking into account \T6\ . we can obtain 

X X ^ x<9a(x,t) ,9a(x,t) *(x,i) ^ 

7Aa(x, i) - *(x, t) + 0/ = U(x, t) ' - —\ XI + V(x, t) ^ - + 9i (74) 

Notice that according to the proposed notations we can rewrite the term {^{x, t) ^"g^f^ — xi in the 
following form: {iJj{x,t)Lf-^a.(x,t) — Lf^4'(x, i)) + (?/'(x, i)LgjQ:(x, t) — Lg-^'i'{x,t))u{x,9,t). Hance it follows 
from (in|) and ^ that ^a{x,t) - *(x,<) + 9i = ip{-)jj)a{x,t). Therefore 9(x,t) = r(-0 + Lp{i;))a(x,t). The 
theorem is proven. 

Proof of Theorem To prove the theorem, first notice that control function H29|) provides the following 
error model dynamics 

^ = -^(V') + z(x,0,t)-z(i,^,i), z(x,0,t) =L^(- e)?A(x,i). (75) 
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By adding and subtracting the function z{i,6,t) from the right-hand side of lf75|l we get the foUowing: 



■ip = - fiip) + z{x, e, t) - z{x, e, t) + z(x, e, t) - z(x, e, t), 

where the difference z{x,6,t) — z{x,6, t) is bounded due to Assumption |51 Denote £{t) = z{x,6,t) — z{x,6,t), 
then 

= -ipiip) + e{t) + z(5c, 0, t) - z{i, e, t), (76) 

where e € L^o- Denote 

f = fi(x)©f^(x)®^f5(x,|,i) 

g = gi(x) ® gj,(x) © 0,.., = (0^__0)^. 

dim X2 

Let us consider the following adaptation algorithm: 

6{x,x,t) = r(0p(x,x,t) r > 0; 0p(x,x,<) = ■i/'(x,t)a(x,t) - ^'(x,t); 

bj = ^(^(x,i))a(x,t) + ^H^_^(x,t)^^^-(V;(x,t)ifa(x,i)-L^vI/(x,i))- 

(-0(x, t)L^a{x, t) - Lg*(x, t)) u{x, hj, 0, t) - A^(i, i), (77) 

where A > 0. Differentiation of fimction Op with respect to time leads to: 

0p(x,x,t) = ?/'(x,i)a(x,t) + V'(x,<)a(x,t) - ^'(x,t) = ?/'(x,i)a(x,<) + V'(x,<)-^"*'^'^^ 9*(x,i) 



9i 



9a(x^ d^{x,t) \ . / ,, 9a(x,^) 9*(M\.' 



9x1 9x1 / \ ' 9x2 9x2 



d\ 9h5 

Taking into account (|77|l and (|23|) we can rewrite H78|l as follows: 

^P(x,x,t) = ^(x,t)a(i,t)-^^^ + 7/^(x,t)^^^ + (7/.(x,t)Lfa(i,i)-if*(ic,t))+ (79) 
(V'(x, i)Lga(x, t) - L~^^{i, t)) u{x, hj, 9, t) + (^^(x, <)^^^ - ^^^) '^'(X' ^) 

Notice also that according to Assumption |S1 

d'^{^,t) ,, 9a(x,0 

and (?/'(x, t) — ■(/'(x, i))Lj^/(x,e)Ct(x, t) G Lqo due to Assumption |51 Denoting (?/'(x, i) — -0(x, i))Lj^/(x,e)Q:(x, t) = 
(5(t) and using equalities ((7^ and lf77|) we can derive that 

e = T{ep + Oi) = r((v3 + (p{^(x, t)))a{x, t) + s{t) - xe), (so) 

where function d{t) is bounded. Let us define the extended state space vector q = x © ^. Furthermore, we 
define 2q(q, 9, t) — z{x, 9, t), aq(q, t) = a{x, t), V-'q(q, t) — ip{x, t). Given the chosen notations, Algorithm (|80|l 
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can be written as follows: 6 = + ip{tpq))cXq{q,t) + S{t) — X9). Moreover instead of equation H76|) we can 

write ■(/; = -ifiip) + Zq{q, 6, t) ~ 2;q(q, 6, t) + £{t). 

It is easy to see that Assumptions O and hold for the extended system. Assumption is also satisfied 
with respect to the goal function due to hypothesis H26|) in Assumption |S1 Indeed, ■!/'q(q, — ip{x,t) e 

Loo X G Loo => ^ G ioo =^ q G Loo- Therefore, according to Assumption |51 and Lemma ^ we can conclude 
that ^(x, i) is bounded and furthermore trajectories x, ^ are bounded as well. Thus property P8) is proven. 
To prove property P9) it is sufficient to notice that 6{t) = either due to the equality da.{x,t) /dx'2 = or 
tp{x,t) = ■0(x, i). Let 9q;(x, t)/9x2 = 0, then P9) follows explicitly from Assumption 1101 and Lemma [5] applied 
to H7t)|l with algorithm 

0(x, i, t) = r(0p(x, i, t) + ei{t)), r > 0; 0p{x, i, t) = V(x, t)a{x, t); 

Bi = ip{ijj{Ti,t))a{Si,t) - ijj{ii,t)da{Si,t)/dt ~ ■il!{x,t)Lfa{i,t) - (i/)(x, t)iga(x, t)) u(x, h^, 0, i). (81) 

which is in fact algorithm (|82|l for A = and \I'(x,i) = 0. If ip{x,t) = ip{x,t), then according to Lemma|21 
algorithm (|77|) with A = 0: 

0(x,x,<) = r(0p(x,x,t) r > 0, 0p(x,x,i) = ^(XjOaCx,*) - *(x,t) 

Bi = .^(V'(x,i))Q:(x,t) + a*(x,t)/5i-V'(x,t)aQ:(x,i)/9t- (-^(x, t)if a(x, i) -Lf1'(x,t)) - 

(V'(x, t)Lga{x, t) - ig1'(i, t)) u(x, hj, 0, t). (82) 

ensures P9) as well. The theorem is proven. 

Proof of Theorem^ We will prove the theorem by induction from order 1 to n for system (|32|l . According 
to the theorem conditions functions /i(xi,0i), ai{xi) are smooth and therefore according to Lemma|21and 
Theorem|21 there exists smooth function u{xi^6i): 

u{xiJi) = -/i(xi,0i) -(pi(V'(a;i)), 0i=7i(0i,p(xi)+0i,/(i)), 71 >0 
Oi.p{xi) = i;{xi)ai{xi) -'^{xi), *(xi) = / V'(a^i) 

~Oi,i = (piiipixi))oti{xi) 

such that tpi{xi,t) S L2 n Loo, ■01 G -^2 for the system of the following type: xi = fi{xi,6i) + u{xi,6i) + £i(t), 
Si(t) G L2. Hence, the basis of induction is proven. 

Let us assume that the theorem statements hold true for the systems of order i, i.e. there exists such smooth 
function Ui{xi, 6i, Vi), Xj, G i?', x^ = (xi, . . . , Xif^ , — (^1, . . . , ^i)^ and the corresponding goal functions 
ipj{xj, t), j = 1, . . . ,i such that -ipji^jjt) G L2 fl Loo, G L2 for system (|32f) of order i: 

Xj ^ fj{xi,...,Xj,ej)+Xj+i, j e {l,...,i-l}, 

Xi = fi{xi,...,Xi,Oi)+Ui+ei{t), ei{t) e L2. (83) 

Therefore, in order to prove the theorem it is enough to show that its statements hold for system H32(l of order 
4 + 1 given that it holds for the systems like (|83() . 
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According to the inductive assumption function Mi(xi, 0^, i/^) is smooth. Then by Hadamar's lemma 
there exists such i^(xj, x-, 0^, 1^, iv^) that Mj(xj, 0^, i^i) - Ui(x^, 0^, 1^, z^j) = i^(xi, x-, 0,,, ^j, iv,j)(xj -x^). Let 
us denote Fj-^j(xi, x^, ^j, i/i) = ||F(xi, x^, 0^, ^j, i/i) |p. Furthermore, due to the theorem conditions functions 
/j(xj, 9j), j = 1, . . . ,i + 1 satisfy the foUowing additional assumptions: 

(/,(x„ 0,) - /,(x;., 0,))' < l|x, - f ^,2(x„ x;.) V e, e ne 

Therefore, it follows from Lemma 13 that there exists a system of differential equations 

^^+1 = f?,+i(€»+i:X,z,iy,+i),| e 

such that 

w,(x,,0j(x,,,|,,0/,j),|,,ivj) - Ui(|,+i,0j(|j+i,|,,0/,j),|,,,i^j) e £2 

Let us introduce new goal function iv^, = Xi+i - 0i(^j+i, 0/,i), i^i) and 

consider its time-derivative ipi+i'- 

i^i+l = /4+l(Xi+l, 0j+l) + - £f^.Mj(^j+i, 0i(^j+i, ^j, 07,j), I,;, I'l) - 

^f«, (1^+1 , , |„ 0/,» ) , |„ ) . (84) 

Denote ei+i{t) = /i+i(xi+i, 6^+1) - fi+l{$^_^_l ®Xi+i,9i+i) and rewrite (|H3Il in the following way: 

^fe^ 0/,,), 1^, I^z + £»+l (0- (85) 

Let us select input u^+i as follows 

Ul+l = -<^,+ l(?A,+ l(x,+ l,|,+ l,|,,IV,;,0/^,;)) +if^.M,(|j+;^,aj(^,^;^,^j,0/^,),|,,Z^j) + 

Lff)U,{^^^^,e,{^,^^,^^,ei.,),t^,i^,) ~ /j+i(^,+i ©a;j+i,a,+i) (86) 

Denoting V'i+i(a;i+i, ^i+i, ^i, I'i, 0/,i) = V'i+i(2;i+i, t) (as are functions of time t) and substitut- 

ing ()86|) and H85|) into 184|) we can write the following expression for ^i+i 

It follows from the theorem conditions that there exits such function Q!i+i(^i+i (Bxi-^-i) that Assumptions |31 and 
Elare satisfied for the function /j+i(^j+i ^i+i)- Consider the following adaptation algorithm: 

Oi+i =7i+i(V'i+i + ¥'i+i(V^i+i(a;i+i,i)))ai+i(|,;+i ®Xi+i), 7^+1 > 
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Realization of algorithms (|88|) is guaranteed by Theorem |21 and can be given as follows: 



dt 

(0) i9a;i4.i 



It follows from Lemma 13 that for the error model H87|l with adaptation algorithm (|88|l and its realization (|89|l 
the following statements hold true: di+i{^i^i(Bx.i+i,t) G iooj "^i+i ^ L2, ipi+i{ipi+i{xi^i,t)) £^2^1 Loo- Given 
that £i+i(t) e -^2 we can conclude that 

Let us denote Ui+i(xj, i>i+i) = Mj+i(xi, iv^+i, iv,,, i>i+i = |j © /^i+i © iv^ © Oj^i. 

According to the introduced notations it is easy to see that statement 2) of the theorem holds. In addition 
to this notice that the choice of appropriate function 1^9^+1 (•) in (|86|) is up to the designer. Therefore choosing 
ipi+i{-) ■■ \ipi+i{-)\ > fcfe+ilVJi+il, > we can guarantee that ?/^i+i(a;j+i, e L2 nioo- 

The last, however, according to the inductive hypothesis implies that ipi^lxi^, t) £ L2 Cl L^o and ip^. G L2 for 
any k — 1, . . . ,i, ^ G L2 C] L^o and ^ G L2. Hence statement 1) of the theorem is proven as well. 

Let us prove statement 3). According to the inductive hypothesis Xi,^j,fi,0i are bounded. Furthermore, 
Oj i{t) is bounded as 0p.i(xi,^j) is smooth function and Oi = "/i{9pi{:iCi,^^) + Bi t)- Then taking into account 
Lemma|31we can conclude that ^^j^i, Vi+i are bounded. Hence I'i+i is bounded. Let us show that x^+i is bounded 
as well. First notice that the difference £i{t) = Wj(xi, 0j(xi, ^i, iv,,) - 1^, ^/^j), /^i) 

is bounded as is smooth and its arguments are bounded. On the other hand we have just shown that 
ipi+i — Xi+i — Ui{(,^^i,9i{$^^^^,(^^,9i^i),$,^,Ui) is bounded. Therefore, Xi+i is bounded. Hence statement 3) is 
proven. 

Derivatives ij^j, j = are bounded as £i{t) is bounded (according to the inductive hypothesis the 

theorem holds for any j = 1 . . . If, however, e{t) is bounded then ipi+i is bounded as well as Ui^i(-), fi+i{-) 
are smooth and x^+i, ^^^j^, f i+i, 0;+! are bounded. Therefore, it follows from Lemma [5] that i/'i+i — > as 
t — > 00. Thus statement 4) is proven. The theorem is proven. 
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